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using Nf = 2 twisted-mass dynamical fermions. The extraction of in the chiral limit re- 
veals the presence of a non-perturbative OPE contribution that is assumed to be dominated by a 
dimension-two (A 2 ) condensate. In this contest a novel method for calibrating the lattice spacing 
in lattice simulations is presented. 
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1. Introduction 

Ajjs is the scale of strong interactions. This parameter has to be taken from experiment and 
can be determined from the running of the QCD coupling constant. This latter can be calculated 
in a variety of non-perturbative ways on the lattice (see ^ [| ||, ||] and references therein). 
In the quenched case ^ the comparison between the perturbative and lattice determinations over 
a large momentum window revealed the presence of a dimension-two (A 2 ) condensate, signaling 
that momenta considered in lattice simulation are in a non-perturbative region. Here we extend the 
strategy of [[7|] to the case of Nf = 2 twisted mass in the sea sector using configurations produced 
by the ETM Collaboration [^], in order to study the effect of the quark mass. 



2. Lattice computation of the coupling in the Taylor scheme 

Following we calculate the strong coupling constant from the ghost-gluon vertex. Gluon 
and ghost propagartors in the Landau gauge are defined as 

G®)^, A) ,( F <»)»V,A) = -^^ (2.D 

where A = aT^ (j8) is the regularisation cut-off. G and F are the gluon and ghost dressing functions 
which can be determined by a non-perturbative renormalization (MOM). In the Taylor scheme [j|], 
where the incoming ghost momentum vanishes, the ghost-gluon vertex does not renormalize. This 
allows for a simple determination of the renormalized coupling constant in this scheme as 

«r(M 2 ) = ^f 1 = lim ^G^.Ay^.A 2 ) ; (2-2) 

in terms of only two-point gluon and ghost dressing function. Here go is the bare strong coupling 
and jit the renormalization scale. This definition can be used in a lattice determination and is to be 
compared with a theoretical formula in order to extract Aqcd- As in the quenched case, using the 
four-loops expression for the coupling constant in the Tayol scheme [R]] 
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where t = In 4V and coefficients are 
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jSo = n--iv> 3 p 1 =p l = m--N f 

p 2 = 3040.48 - 625.387 N f + 19.3833 Nf 
j8 3 = 100541 - 24423.3 N f + 1625.4 Nf - 27 .493 Nj , (2.4) 

Extracting A^ from the lattice data at each /I 2 using this perturbative formula does not lead to 
a constant value. To understand the mismatch beetween lattice and perturbative determination, a 
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non-perturbative OPE correction to the perturbative formula is to be considered. This accounts for 
the minimal power correction associated to the presence of a dimension-two (A ) condensate: 

where Aq C d is some perturbative scale. This will cure the mismatch and lead to a good 
determination for At in the Taylor scheme, which eventually can be be related to the value of the 
scale in the MS scheme through 

ci _ 507 -4(W/ 

^ms = 'lb = ~192-4SN f = . 541449 . (2 6) 

At 

3. Artefacts 

We exploited data from ETMC configurations at maximal twist for a variety of run parameters 
(tab. |]) in order to study physical and systematic effects in our determinations. This have the 
main advantage of reducing the discretization artefacts to 0{a 2 ), where a is the lattice spacing. 
Nevertheless, artefacts are expected to came at different levels. A first kind of artefacts that can 



be systematically cured Jll] , |12| are those due to the breaking of the rotational symmetry of the 
euclidean space-time when using an hypercubic lattice, where this symmetry is restricted to the 
discrete H(4) isometry group. It is convenient to compute first the average of any dimensionless 
lattice quantity Q(ap^) over every orbit of the group H(4). In general several orbits of H(4) 
correspond to one value of p 2 . Defining the H{4) invariants = Ym=i Pu' ^ me lattice spacing is 
small enough such that e = a 2 p^/p 2 « 1, the dimensionless lattice correlation function can be 
expanded in powers of £: 

Q(a 2 p 2 y p W,a 6 p^,a 2 A 2 QCD ) = Q(a 2 p 2 ,a 2 A 2 QCD ) + ^ " ^ 



a 2 ^ T + --- (3.1) 
o P 2 



H(4) methods are based on the appearance of a &{a 2 ) corrections driven by a p^ term. The basic 
method is to fit from the whole set of orbits sharing the same p 2 the coefficient dQ/de and get the 
extrapolated value of Q, free from H{4) artefacts. 

A second kinf of artefact could come from dynamical quark masses. We will argue that this is 
a ff{a 2 ix 2 ) effect and therefore that it is a lattice artefact. We have calculated the //(4)-free ghost 
and gluon dressing functions and combined in order to calculate the //(4)-free lattice coupling 
through eq. (2.2). In Fig. || one can see the Taylor coupling after hypercubic extrapolation for 



different pL q at fixed j8 = 3.9 and 4.05. Indeed, a dependence in \i q is clearly seen. If it is an artefact 
the dependence should be in a 2 jl 2 . If it is an effect in the continuum it should be some unknown 
function of the physical mass \i q . Trying an &{a 2 pL 2 ) dependence, we write the expansion : 



a T (a 2 P 2 ,a 2 n 2 ) = a T {p 2 ) +R (a 2 P 2 ) a 2 ]x 2 R (a 2 p 2 ) = J^ 1 \ (3.2) 



Provided that the first-order expansion in eq. ( |3.2[ ) is reliable, a linear behaviour on a \l q has to 



be expected for the lattice estimates of aj for any fixed lattice momentum computed from simula- 
tions at any given j8 and several values of \L q . We explicitely check this linear behaviour to occur 
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for the results from our /3 = 4.05 and j3 = 3.9 simulations and show in Fig. |3| some plots of Sj- 
computed at j8 = 4.05 (where four different quark masses are available) for some representatives 
lattice momenta in terms of cP-\i 2 . In fig. |J we plot Ro(a 2 p 2 ) as a function of ap computed for the 
four lattices simulations at j8 = 4.05 with different quark masses and for the three ones at /3 = 3.9. 
Indeed, it can be seen that a constant behaviour appears to be achieved for p > p m [ n ~ 2.8 GeV. We 
will not risk an interpretation of the data below (ap) m i n . The striking observation here is that above 
p m m both lattice spacings exhibit a fairly constant Ro(a 2 p ) and a good enough scaling between 
both j3's. The fact that Rq with our present data goes to the same constant for both /3's, leads us to 
consider that the /J. dependence of a is mainly a lattice artefact (else it should be a function of jJ. 
and not of apL). 

The main result of this work is taking into account the effects due to dynamical quarks in a 
global analysis of the lattice determinations. This lead to a proper extrapolation to the continuum 
limit, which can be compared with continuous formula in order to extract A^j§-. 
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Figure 1: Left: Run parameters of the exploited data from ETMC collaboration. Right: Best-fit parameters 
for the ratios of lattice spacings, A^g and the gluon condensate (for which a(3.9)qo = 4.5 is chosen). For 
the sake of comparison, we also quote the results from [ jl3| ] that were obtained by computing the hadronic 
quantity, ro/a(j3), and applying to it a chiral extrapolation. 



4. Ajyfg and the gluon condensate 



The running of (Xt given by the combination of Green functions in eq. ( |2.2| ) and the extrapo 



lation through eq. (p.2|), provided that we are not far from the continuum limit and discretization 
errors are treated properly, depend only on the momentum (except, maybe, finite volume errors at 
low momenta). The supposed scaling of the Taylor coupling implies for the three curves plotted 
in fig. f| to match to each other after the appropriate conversion of the momentum (in x-axis) from 
lattice to physical units, with the multiplication by the lattice spacing at each j3. Thus, we can apply 
the "plateau"-method described in [||] for the three /3 's all at once by requiring the minimisation of 
the total x 2 : 

2 



a(j3i) a(j8 2 



*i^«ffl>m = &* — % — - (41) 
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p=3.9 




Figure 2: Left: The Taylor couplings estimates, after 7/(4)-extrapolation, at j3 = 3.9 for fi q = 
0.004,0.0064,0.010. Right: The slopes for the mass squared extrapolation in terms of ap computed for 
the four lattices simulations at j8 = 4.05 (32 3 x 64) with a\l q = 0. 003,0.006,0.008,0.012 and for the three 
ones at j3 = 3.9 (24 3 x 48) with ajx q = 0.004,0.0064,0.010. 




Figure 3: We plot the values of the Taylor coupling at j3 = 4.05, computed for some representative values 
of the lattice momentum, a(4.05)p = 1.08, 1.18, 1.24, 1.36, 1.45, 1.52, in terms of a 2 (4.05)jl 2 and show the 
suggested linear extrapolation at a 2 j^ q = 0. 



where the sum over j covers the sets of coupling estimates for the three /3's (j8o = 3.9, j8i = 4.05, 
/?2 = 4.2), the index i runs to cover the fitting window of momenta to be contained in a region in 
which the slope Rq ~ —90 was found to be constant. A,-(j3y) is obtained for any j8 ; by requiring the 
best-fit to a constant; c results from the best-fit: it is the Wilson coefficient of the gluon condensate 
in eq. (2.5), where the leading logarithm correction is now taken into account, where a(/3o)<7o = 4.5 
(this means q$ 10 GeV) was chosen. The function % 2 is minimised over the functional space 
defined by the four parameters that are explicitly put in arguments for eq. ( Ll)'s l.h.s.: a(po)A-^, c, 
alfo)' o(Bo) " Thus we obtain all at once A^- and the gluon condensate, in units of the lattice spacing 
for j3o = 3.9, and the ratios of lattice spacings for our three simulations after the extrapolation to 
the limit \i q — > (see tab. [[]). The errors are calculated again by jackknife analysis. The ratios of 
lattice spacings can be applied to express the momenta for all the three sets of coupling estimates 
plotted in fig. ^ (left) in units of the lattice spacing at /3 = 3.9. Thus they indeed match each 
other and fit pretty well to the analytical prediction with the best-fit parameters for A^ and the 
gluon condensate, in units of l/a(3.9) (see tab. [j]), as can be seen in the plot of fig. ^| . A detailed 
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a(p)p ' a(p=3.9)p 



Figure 4: Left: The Taylor coupling, free of H(4) and mass-quarks artefacts, for the three j3 = 3.9,4.05,4.2 
and plotted in terms of the lattice momentum a(j3)p. Right: The scaling of the Taylor coupling computed 
by for the three j3 = 3.9,4.05,4.2 is shown. The lattice momentum, a(j5)p in the x-axis, is converted to a 
physical momentum in units (the same for the three j3's) of a(3.9) _1 . 



discussion about systematics can be found in [jl4J indicating that main sources of errors are under 



control. Assuming the value a(3.9) = 0.0801(14) fm [13], we quote our result as 
^MS ~ (330 ± 23) x MeV , ^)<A% = (2.4 ± 0.8) x (^"^ — 2 



5. Conclusions and outlooks 

We computed the renormalized strong coupling constant analyzing a variety of Nf = 2 gauge 
configurations generated in the ETM Collaboration. We performed an elaborated treatement of the 
lattice artefacts and a precise estimate of the couplings at the infinite cut-off limit. The coupling 
estimates for lattices at different j3's were seen to match pretty well, as should happen if the cut-off 
limit is properly taken, when plotted in terms of the renormalization momenta converted to the 
same units by applying the appropriate lattice spacings ratios. These ratios could be either taken 
from independent computations or obtained by requiring the best matching with pretty compatible 
results. Thus, once we are left with the estimates of the coupling constant extrapolated at vanishing 
dynamical mass \L q , for every value of the renormalization momentum, jx, they were converted via 
a fit with a four loops formula into the value of A^-. As in th quenched case a condensate (A 2 ) 
is needed in order to get a constant A^. As an outlook, we want to apply the same analysis to 
the case of lattice QCD with Nf = 2 + 1 + 1 and N + / = 4 dynamical flavors. This will lead to 
give a reliable lattice prediction for the coupling constant, say at Mz, to be compared with available 
experimental determinations. 
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